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Abstract 

We present a phenomenological model of melting in nanoparticles with facets that are only 

C^ ' partially wet by their liquid phase. We show that in this model, as the solid nanoparticle seeks 

to avoid coexistence with the liquid, the micro canonical melting temperature can exceed the bulk 

C I melting point, and that the onset of coexistence is a first-order transition. We show that these 

■ results are consistent with molecular dynamics simulations of aluminum nanoparticles which remain 

CN . solid above the bulk melting temperature. 
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The thermodynamics of small systems is becoming increasingly relevant technologically 
as the size of devices enters the nanoscale [l[. A classic example of a phase transition in 
a small system is the melting of a free nanoparticle. Generally, one expects the melting 
temperature to decrease in proportion to the surface area to volume ratio of a particle 
because the surface free energy of a molten droplet is less than that of the corresponding 
solid particle. This scaling has now been confirmed in many instances J2; S |J|, although 
a number of exceptions to this general rule have appeared. For instance, as the sensitivity 
of free nanoparticle calorimetry has improved [5|, a so-called 'non-scaling' regime has been 
observed, where the melting temperature is seen to varverratically with particle size p|. 
Further, small clusters of both lead |7| and gallium 8|, |9| have been observed to melt at 
temperatures higher than the bulk melting temperature, T^. 

Superheating has also been seen in the surface melting (SM) of bulk crystals [lO|. Surface 
melting is generally thought to occur on metal surfaces for which A7 = 7^^ — 7^/ — 7zt, > 
(where 7^^ is the solid-vapor (SV) interfacial energy, 7^/ is the sohd-hquid (SL) interfacial 
energy and 7;^, is the liquid-vapor (LV) interfacial energy) i.e. below Tc a liquid layer on 
the surface is frequently observed 11|, the thickness of which diverges logarithmically as 
T — >■ Tc. However, for particularly low-energy metal surfaces with A7 < (such as Pb (111) 



]J, Al (111) [l^ and (100) [JJ, and Au (111) [15|) surface melting is not seen, and often 
the surface can be heated above Tc |lo|. Such surfaces are called non- melting (NM). Indeed, 
this is thought to be the reason why supported Pb nanoparticles with only NM (lll)-facets 
can be superheated several degrees above the Tc |2|, although in this case the superheated 
solid cluster is presumably only metastable |lfil |. 

In an isolated nanoparticle where the internal relaxation time is less than the time for 



equilibration with the environment 
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, the situation is somewhat different. Such a particle 



will follow a microcanonical caloric curve [13| , where the transition between solid and liquid 
phases will generally occur via solid-liquid phase coexistence, even if the particle only has 
NM facets. In a small particle, the cost of forming the solid-liquid interface during phase 
coexistence is comparable to its total energy, and this can lead to an S-bend in the micro- 
canonical caloric curve |19|, negative heat capacities p, the avoidance of phase coexistence 
entirely |20] or dynamic coexistence between fully solid and fully liquid states [21]. Phase 
coexistence in clusters with NM facets will be particularly unfavorable and, while coexis- 



tence has been observed in molecular dynamics (MD) simulations of such clusters 22], there 



is evidence to suggest that the transition between the sohd and sohd-hquid coexisting state 
in this case is first-order rather than continuous J20|, |23| . 

In this letter, we use microcanonical critical droplet (MCD) theory |24] to consider melting 
and solid-liquid phase coexistence in free clusters with NM surfaces. Interestingly, we find 
that not only can superheating of such solid clusters occur, but in this case, the superheated 
clusters are stable rather than metastable. We also find that the transition from solid to 
solid-liquid coexistence is first-order for NM surfaces in MCD theory, and that there is a 
minimum size below which solid-liquid phase coexistence is unstable. Both these findings 



are in agreement with results from previous MD simulations J2Q|]. Finally, we conduct new 
MD simulations of solid-liquid coexistence in Al truncated octahedra (TO) with NM (111) 
and (100) facets, where we observe both superheating and a first-order transition to solid- 
liquid coexistence in agreement with the theory. In the MD simulations we use a glue 
potential developed for Al by the first-principles force- matching method J25|. In particular 
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so we 



this potential has been used to study melting of the Al (111) and (100) surfaces 
expect it to be applicable here. 

For a coexisting cluster of radius R with NM surface facets, the melt will not fully wet 
the solid so the usual assumption of a spherically-symmetric liquid layer (see Fig. Q) is 



inappropriate. One approach is to use a geometry of two overlapping spheres J22] but such 
a model requires a two parameter description. Here we will use a simple one-parameter 
model geometry to describe the NM case (Fig^. Specifically we characterize the coexisting 



cluster by the height /i of a spherical 
geometry seen NM Pb nanoparticles 



iquid cap (we note this does resemble the solid-liquid 
20]). In this model the cluster is fully solid at /;, = 



and fully liquid at /i = 2R. The total surface energy of the cluster F may now be written as 

F = TT (2i? {2R - h) 7,, + Rh^i, + h{2R- h) 7,,) . (1) 

Note that we have neglected any interaction between the LV and SL interface. For short 
range interactions this interaction decays exponentially (ll| with a correlation length ^, so 



F wi 

a|i^|. 



1 be well approximated by (^ for /i > ^. For Al (111), ^ has been estimated to be ~ 3 



Now consider the dependence of F on h. It is easy to show that —^ < only when 

h>h* = -—R. (2) 
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Thus for a cluster with NM surfaces (A7 < 0), the total surface energy can only be reduced 
by increasing the liquid volume for a liquid cap thickness h > h* > 0. This is in contrast to 
the SM case, where F is always reduced by increasing the liquid volume. The inequality in 
((21) strongly suggests that any coexisting state in a cluster with NM surfaces will be unstable 
for h < h*. We note that with A7 < 0, in the spherically symmetric case (Fig. ^ one can 
show that the surface energy must increase with any increase in liquid fraction, so there is 
no critical liquid fraction. Nonetheless, a critical condition such as Q does not appear to 
be a peculiar feature of the spherical cap geometry. For instance, if one considers a solid 
sphere in contact at a point with a liquid droplet, one can again show there is a critical liquid 
fraction where surface energy will decrease with an increase in liquid fraction. Indeed the 
existence of a critical liquid fraction is consistent with the first-order nature of the transition 
to coexistence seen in simulations of Pb clusters with NM surfaces J20|. 

To make the discussion more concrete, we construct a caloric curve for a cluster with NM 
surfaces us,,, MCD tKeo.y fl Q. TKe total ene,, E of t>,e cluste. ,s t>,e su. of t>,e 



volumetric and surface contributions: 

E={V-Viih))es + Viih)ei + T{h) (3) 

where V = 47ri?^/3, Vi{h) = n/3h'^{3R — h) (assuming that the density is unchanged at 
melting) and es{i) is the energy density of the solid (liquid). The total entropy, Sm, at E is 

S^{E, h) = {V- Vi{h)) ss{es) + Vi{h)si{ei) (4) 

where Ss{es) = s^ + clog f ^'^^ + 1 j is the entropy density of the solid region, and similarly 
si{ei) is the entropy density of the liquid. Here c is the heat capacity (assumed to be the 
same in the solid and liquid states) and quantities with a superscript c indicate values at 
Tc- At a given total energy _E, the state of the coexisting cluster is described by the value of 
h = h'^ that extremizes the entropy Sm- The onset of coexistence occurs at an energy when 
Smih^E) = Ss{E) and complete melting occurs when Sm{h),E) = Si{E). These equations 
must be solved numerically to find the energies at which the transitions occur. 

In Fig. 121 we plot the caloric curve constructed by extremizing Q for a 3.5 nm Al cluster 
(which corresponds to 9590 atoms). We have assumed the cluster has only (111) facets 



and have used material parameters estimated using the glue potential for Al |l 



particular these parameters give A7 = -2.3 meV/A^ and -^ ~ —0.23 [l6[, so in a 3.5 nm 



Q, 



25|. In 



cluster, h* ~ 8A > ^. Thus, unless the interaction between the LV and SL interfaces causes 
A7 to change sign, our discussion will be unaffected by the approximation in (Q). Note 
the first order character of the transition at the onset of coexistence, consistent with the 
inequality Q- In this case we calculate the the onset of coexistence occurs at h=19.8 A > h*. 

We also note the superheating of the solid evident in Fig. |21 the cluster does not start 
to melt until T > Tc. In the SM case, where A7 > 0, a solid cluster begins to melt below 
Tc because it can convert an increment of surface energy into latent heat. In the NM case, 
the cluster will not melt until it is favorable to melt a finite volume with h > h*. With the 
parameter values used here, the temperature of the solid at melting satisfies T > Tc. Indeed, 
it is important to emphasise that the superheated state here is stable rather than metastable. 
We note however, that for smaller values of A7 or L, for instance, we have found that the 
melting temperature can drop below Tc. Thus superheating is not a necessary consequence 
of the model but results from the particular parameters used for the Al (111) surface. 

With these parameters, for R < 3.4 nm we find that the is no energy E where Sm{h'^ ■, E) > 
Ss{E) and Sm{h'^,E) > Si{E). In other words the coexisting state is always unstable and 
the cluster will melt fully without undergoing coexistence when Ss{E) = Si{E). At this 
point it is useful to note that the model geometry in Fig. Q probably overestimates F for the 
coexisting particle. Thus we expect the theory here will underestimate the stability of the 
coexisting state, and overestimate the degree of superheating. In addition, it is important 
to point out that our model neglects the size-dependence of properties such as the latent 
heat of melting, L (which appears implicitly in (jlj) as Sf — Sg = L/Tc). 

However we now give an example of this superheating and the corresponding first-order 
transition to coexistence obtained by MD simulations of Al nanoparticles using the glue 
potential J25|. With this potential, closed-shell Al clusters with 586 or more atoms favor 
TO structures 26|. We simulated caloric curves for the 586, 1289, 2406 and 4033-atom TO 
clusters. At each energy, the cluster was equilibrated for 0.34 ns (with time step 1.35 fs), 
then the kinetic energy was averaged over a further 0.34 ns. The energy increment between 
simulations was 1.0 meV/atom with energies adjusted by a uniform scaling of the kinetic 
energy. Coexisting states were identified by the appearance of a bimodal distribution of 
atomic mobilities J22,|23[. 

The caloric curve for the 4033 atom cluster is shown in Fig. El Firstly, we note that the 
transition from the solid to solid-liquid coexistence is first-order, as we expect from (J2j, with 



a large drop in temperature at the onset of coexistence. Further, the cluster melting begins 
at 975 K which considerably exceeds the bulk melting temperature {Tc = 939 K using the 
glue potential). Thus, qualitatively, the 4033-atom cluster shows similar behavior to that 
predicted by MCD theory. 

Of course, the overheated solid could just be metastable rather than stable as suggested 
by MCD theory. To test this, we selected a cluster state at a point in the coexistence region 
and removed energy at the corresponding rate. The hysteresis in the cluster state is evident 
in Fig. El but the cluster fully solidifies at a temperature of 944 K. This shows that the 
superheated solid is stable relative to the coexisting state over a range of energies. Thus 
we conclude that the melting temperature of the 4033 cluster lies between 975 and 944 K. 
Interestingly, the 4033-atom cluster displays two distinct coexisting states shown in Fig. |^ 
In the first state, which is at least metastable, melting is confined largely to two neighboring 
(lOO)-facets and the (lll)-facet adjoining them. In the second state, all but one of the 
(lOO)-facets have melted. 

Fig. El compares the melting points predicted by MCD theory with those obtained by 
MD simulation. Also indicated on the plot is the threshold for coexistence predicted by the 
theory. In the MD simulations only the 2406 and 4033-atom clusters exhibited static phase 
coexistence prior to melting. No indications of static, transient or dynamic coexistence were 
seen prior to melting in the simulations of the 586- atom cluster although the 1289- atom 
cluster displayed transient coexistence (i.e. dynamic coexistence between the solid and a 
solid-liquid state, suggesting that the coexisting state is metastable). Thus the threshold 
for the stability of coexistence lies between the 1289-atom and the 2406-atom clusters. This 
is a smaller threshold than we obtained from MCD theory, but as noted earlier, our model 
geometry (Fig. P) is likely underestimate the stability of the coexisting cluster. We note that 
MCD theory predicts that the superheating effect will peak at a cluster radius oi R = 5.5 
nm. As -R — > oo, the melting temperature approaches Tc. 

In previous work, we have constructed microcanonical caloric curves for Ag, Cu, Ni (23l |. 
^d J22| and Pb j2fl] nanoparticles using molecular dynamics with embedded atom method 
28[ potentials. Our initial studies of Pb icosahedra with NM (lll)-facets 12] found that 
in particles with radii less than 2 nm, solid-liquid coexistence was absent prior to melting, 
and that at sizes above this the transition to coexistence was first-order. We did not see 
superheating in the Pb icosahedra - we recall this is not a necessary consequence of MCD 



theory. In Ag, Cu, Ni and Pd we saw coexistence in much smaller particles. Furthermore, the 
transition to coexistence in these clusters was always continuous. In Ag we saw evidence for 
a coexisting solid-liquid state in a 309-atom cluster, realized as dynamic coexistence between 
a solid-liquid state and a fully solid state. Similarly in Pd, we saw a coexisting solid-liquid 
state in a 309-atom cluster, although it was realized as dynamic coexistence between a solid- 
liquid state and a fully liquid state. Static solid-liquid coexistence was observed in both Ni 
and Cu 561-atom icosahedra. These previous findings are consistent with our considerations 
here, where clusters with SM facets exhibit a continuous transition to coexistence down to 
sizes of -R ~ 1 nm, and clusters with NM facets avoid coexistence in sizes below R ^ 2 and 
exhibit a first-order transition to coexistence at sizes above this. Thus the stability of phase 
coexistence in a cluster is strongly dependent on the sign of A7. 

In conclusion, we have found theoretical evidence for superheating in clusters with NM 
surface facets in the microcanonical ensemble. The superheating is associated with a mini- 
mum stable liquid volume, which also requires a first-order transition at the onset of melting. 

n 

Unlike previous observations of a metastable superheated state supported clusters [3], the 
microcanonical superheated state discussed here in free clusters is expected to be stable. 
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FIG. 1: Model geometries for coexisting clusters. For clusters with SM facets (A7 > 0), the 
spherically symmetric model on the left can be used. For clusters with NM facets (A7 < 0), we 
use the spherical cap model on the right. 
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FIG. 2: Caloric curve for a 3.5 nm Al cluster constructed using MCD theory. 
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FIG. 3: Caloric curve for a 4033-atom Al cluster. The curve was constructed both by heating 
from the solid state (squares) and cooling from a coexisting structure (circles) that emerges as the 
cluster approaches the liquid state. 





FIG. 4: Snap shots showing the state of the 4033-atoni cluster at E = —2.966 eV/atom (left) and 
E = —2.960 eV/atom (right) during the cooling phase of the caloric curve in Fig. |31 Light grey 
atoms are of low mobility (solid) and the dark grey atoms are of hight mobility (liquid). 
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FIG. 5: The dependence of melting point on size as calculated using MCD theory for an Al 
(lll)-faceted particle (solid lines), and several data points calculated from MD simulations. 
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